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Neutral dynamics, where taxa are assumed to be demographically equivalent and their abundance 
is governed solely by the stochasticity of the underlying birth-death process, has proved itself as 
an important minimal model that accounts for many empirical datasets in genetics and ecology. 

However, the restriction of the model to demographic [O(VJl)] noise yields relatively slow dynamics 
that appears to be in conflict with both short-term and long-term characteristics of the observed 
systems. Here we analyze two of these problems - age size relationships and species extinction time 
- in the framework of a neutral theory with both demographic and environmental stochasticity. 

It turns out that environmentally induced variations of the demographic rates control the long¬ 
term dynamics and modify dramatically the predictions of the neutral theory with demographic 
noise only, yielding much better agreement with empirical data. We consider two prototypes of 
’’zero mean” environmental noise, one which is balanced with regard to the arithmetic abundance, 
another balanced in the logarithmic (fitness) space, study their species lifetime statistics and discuss 
their relevance to realistic models of community dynamics. 

PACS numbers: 87.10.Mn,87.23.Cc,64.60.Ht,05.40.Ca 


INTRODUCTION 

The theory of population and community dynamics 
is a central mathematical tool in many branches of life 
science, including genetics, ecology and evolution. This 
broad field of research is dominated by two competing 
and complementary paradigms. Darwinian natural selec¬ 
tion suggests that the fitness of phenotypes and species 
is different and stresses deterministic effects like the sur¬ 
vival of the fittest, downplaying the role of noise and 
fluctuations. On the other hand, neutral theories as¬ 
sume that selective effects are relatively weak and dif¬ 
ferent taxa/haplotypes admit almost identical htness, so 
that the main driver of population dynamics is stochas¬ 
ticity. 

Within the neutral framework, hrst Kimura’s theory 
of molecular evolution [1] and more recently Hubbell’s 
universal neutral theory of biodiversity (UNTB) [5H3] 
have both attracted a lot of attention. In the latter case, 
the successful explanation of empirical species abundance 
distribution curves by a simple theory with only two pa¬ 
rameters [3i appears as a very appealing minimalis- 
tic model, especially when compared to niche-based ap¬ 
proaches that usually require the reconstruction of many 
parameters (such as the relative fitness of species), a very 
difficult task in high-diversity assemblages [7]. 

These neutral models assume that the main driver of 
community dynamics is demographic stochasticity, i.e., 
the noise embodied in the birth-death process of individ¬ 
ual agents, with (if the whole community has to keep a 
fixed size) the expected number of descendants for each 
individual being precisely one. Accordingly, a population 
of N individuals will produce N offspring on average, 
and the per-generation fluctuations will be proportional 


to -y/fV, which is quite a weak noise in the limit of large 
N. This restriction of the neutral model to pure demo¬ 
graphic noise leads to a few severe difficulties when its 
predictions are compared with empirical patterns. 

On evolutionary time scales, the two main unsolved 
problems are the age-size relationships and species ex¬ 
tinction time. In a neutral theory with pure demographic 
noise both the age (measured in generations) of a species 
and its time to extinction are proportional to its abun¬ 
dance. This timescale is ridiculously long for, inter alia, 
various species of trees (with generation time of 50y) 
and for passerine birds (generation time 3y), as noted by 
many authors [HHII]. The fossil data, which indicate that 
species lifetime is typically a few million years, is again 
in contradiction to the N generation estimate for com¬ 
mon species. As Robert Richiefs summarized his find¬ 
ings [3], “drift is simply too slow to account for the rate 
of turnover of passerine birds”. The idea of protracted 
speciation [T^] , suggested to account for the apparent un¬ 
derrepresentation of rare species and their relatively long 
lifetime, cannot resolve these difficulties. 

On the ecological time scale, empirically observed fluc¬ 
tuations in abundance are usually too strong to be ex¬ 
plained by UNTB [T5HTS] . Moreover, UNTB cannot ex¬ 
plain the scaling of fluctuations variance with population 
size: with pure demographic noise the theory predicts a 
linear scaling but empirical analyses show a prevalence 
of super-linear dependence m The decay of commu¬ 
nity compositional similarity is again much faster than 
the UNTB predictions [T5] . 

A very plausible (and one in any case necessitated by 
biological reality) generalization of the neutral theory 
that may resolve many of these problems while preserving 
the minimalistic character of the model, is to add envi- 
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ronmental stochasticity to the dynamics. In nature the 
demographic response of individuals to environmental 
variations is (at least partially) correlated within species, 
while different species may show transient fitness advan¬ 
tage at various times due to differences in their temporal 
niches. In our proposed generalized neutral model all in¬ 
dividuals are demographically equal on average but the 
relative fitness of a population fluctuates in time. Envi¬ 
ronmental stochasticity generates 0{N) short-term fluc¬ 
tuations in population abundance, closer in size to those 
observed in reality. 

To provide an order of magnitude estimate for the 
strength of environmental noise in empirical systems, the 
variance through time plot of the logarithmic abundance 
ratio is presented in Fig. [^for three empirical datasets, 
showing the linear increase characteristic of environmen¬ 
tal stochasticity. As discussed in the slope of this 
curve is proportional to the strength of environmental 
variations. In the discussion section we plug this number 
into our mathematical expressions to show that realistic 
noise indeed solves the UNTB timescale problem. 

Our hrst major result, then, is that a neutral model 
with both demographic and environmental noise may 
solve the mismatch between the predictions of the purely 
demographic UNTB and the empirical evidence, short¬ 
ening significantly the number of generations needed for 
the ancestry of a single founder to reach large abundance. 
The linear scaling of time (in generations) with the abun¬ 
dance N is replaced by logarithmic scaling, yielding rea¬ 
sonable age-size estimations. Our second result is that 
the time to extinction of an abundant species is short¬ 
ened significantly. As described below, the various possi¬ 
bilities obtained for species lifetime statistics correspond 
to previously, evidence-based, suggestions. 

The implementation of environmental stochasticity 
into a neutral model poses an important conceptual prob¬ 
lem. As noticed a while ago m there are two differ¬ 
ent scenarios for balanced (zero mean) environmental 
stochasticity: the balance may be either arithmetic (i.e., 
the noise statistics is such that the mean abundance is 
kept fixed; for example if during a bad year the popu¬ 
lation shrinks by 1/2, during a good year it grows by 
50%) or logarithmic (in a bad year it shrinks to 1/2, say, 
and in a good year the population doubles). These two 
scenarios are analogous to the use of Ito or Stratonovich 
calculus in the white noise limit [33l [34] . 

For a fixed-size community, where different species are 
playing, more or less, a zero sum game, none of these sce¬ 
narios provide a satisfactory description of the dynamics. 
Under arithmetically balanced noise all species species 
eventually go extinct with probability one, while under 
log-balanced noise the size of the community grows ex¬ 
ponentially, both in contrast with the fixed size require¬ 
ment. Introducing environmental stochasticity into the 
UNTB is a delicate task, in which the interplay between 
environmental and demographic fluctuations and other 


effects should be taken into account. Here we do not 
solve this problem, although we provide some prelimi¬ 
nary considerations in the discussion section. 

However, any model with balanced noise must be 
somewhere between the log-balanced and the arithmetic- 
balanced extremes considered here and, as will be shown 
below, in both cases the abundance-age relationships are 
logarithmic so the species lifetime problem is solved in 
any case. On the other hand, the log and the arithmetic 
dynamics differ dramatically with respect to species life¬ 
time: the tail of lifetime statistics predicted by the 
purely demographic theory is replaced by a slower 
(^-3/2) jjecay if the noise is log balanced and by much 
faster, exponential decay, for arithmetically balanced 
stochasticity. In the literature one may hnd empirical ev¬ 
idence for all these behaviors [231121| ; the combination of 
demographic and environmental stochasticity is endowed 
with the required flexibility to account, in the appropri¬ 
ate regime of parameters, for many observed patterns of 
species/genus lifetime statistics. 


THEORETICAL FRAMEWORK 

Demographic noise reflects variability in reproductive 
success which is uncorrelated among individuals in the 
population. This noise is fully characterized by 
the probability of an individual to produce n offspring 
during its lifetime. The average number of offspring is 
~ ^o, and for a fixed size population Ro = 1. 
Another important parameter is the variance of the num¬ 
ber of offspring tr^, which measures the differences in 
reproductive success among individuals: in animal and 
human populations, where the number of offspring is be¬ 
tween zero and 10, say, and decays sharply with n, 
the variance is of order unity; for a virus population, 
where some viruses infect a cell and produce 10000 off¬ 
spring while the others die childless, the variance is much 
larger. 

The theory of a fixed-size population under pure de¬ 
mographic noise traces back to the work of Galton and 
Watson [5S|. In Appendix A we present a generating 
function analysis of this case. Species lifetime statistics 
are reflected in 5(t), the chance of the lineage of a single 
individual to survive until t (time is measured in units of 
generations), which is shown to satisfy 


This implies that the survival probability of a species 
decays at long times like \/t and that the statistics of 
species lifetime admits a 1/t^ tail [33]. Since the average 
abundance of a species is constant, the fact that only a 
fraction of the species survive by time t implies 
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m=0.00414181 Py^i=0.0166964 R^=0.966886 m=0.0549023 P^^i=<0.00001 R^=0.722689 m=0.0120791 P^gi<0.00001 R^=0.972936 



At [5 years] At [years] [years] 


FIG. 1: Variance through time plots for three empirical datasets: A) Trees of the tropical forest in Barro Colorado Island 
(BCI), considering trees with diameter > 10mm (data are from www.ctfs.si.edu) from 1990 to 2005 (every 5 years). B) Fish 
samples collected from the cooling -water filter screens at Hinkley Point B Power Station in Bristol, from 1980 to 2008 m- 
We only consider estuarine fish species (not crustaceous organisms). C) herbaceous plant dataset [50] comprising a series of 51 
quadrats of 1 m^ from mixed Kansas grass prairies where all individual plants were mapped every year from 1932 to 1972. In 
all cases we gathered log{Nt+At/Nt) from every two same species abundance observations that are At apart, and plotted the 
variance of this quantity versus At (the details of this technique are presented in [15]). m denotes the angular coefficient of the 
linear fit of data points, while Pvai and give the significance and the goodness of the fit, respectively. In all three cases the 
variance appears to grow linearly in time, as suggested for a system with logarithmically balanced environmental noise, where 
the slope m indicates the value of the effective D defined in The intercept of the linear fit is above zero for the fish and 
the grassland, indicating the effect of sampling errors m The negative intercept for trees indicates a delayed response of the 
system to the changing environment. Since the timescales are relatively small, the same analysis for arithmetically balanced 
noise, where the variance of {Nt+At — Nt)/Nt is plotted against t, yields almost the same results (not shown). 


that their mean abundance should grow linearly with t, 

2 

Nsurvit) - yt. (2) 

where Ngurvit) is the average population abundance at 
t, conditioned on survival. These results demonstrate the 
two main problems of a neutral theory with pure demo¬ 
graphic noise, when confronted with empirical findings 
and a priori considerations: first, the time needed for 
a species to reach abundance N (starting with one in¬ 
dividual, i.e., point speciation) is typically of order N 
generations. Second, the typical time to extinction of 
a species of abundance N is again of order N genera¬ 
tions, since the theory is neutral and the ancestry of any 
individual evolves independently (Formally the average 
time to extinction diverges due to the 1/t^ tail). Slight 
modification of the model, like implementing a zero-sum 
game in the community (such that the size of the com¬ 
munity is kept fixed in the strong sense, as opposed to 
keeping the average size fixed by assuming Ro = 1 ), do 
not significantly change these conclusions m- 

Adding environmental stochastieity to this model im¬ 
plies that, as the environmental conditions vary, the de¬ 
mographic success of the whole population varies accord¬ 
ingly, so the average reproductive success, i?o, becomes 
time dependent. In some cases, discussed in [18], the ad¬ 
dition of environmental noise leads to a stabilization of 
the populations around some equilibrium value because 
of the (quite counterintuitive) storage effect. This paper 


deals with models that have no such effect, a set which 
includes any environmental noise generalization of the 
dynamics considered in [^. 

Environmental conditions admit some correlation time 
T. Speaking about a good/bad year (in terms of precip¬ 
itation, winds etc.) one assumes that the environment 
(in the general sense, including the effects of competition 
with other species) was, in general, favorable or hostile to 
a specific species during this period. If the demographic 
rates are kept fixed during T and the abundance of a cer¬ 
tain species at time s is A^s, then typically Ng+i = e'^^Ng 
where 7 ^ = (i?o — 1)T. The simplest way to define a ’’bal¬ 
anced” environmental noise is to assume that 7 ^ (the fit¬ 
ness parameter, or the deviation of Ro from unity) is an 
identically distributed random variable with zero mean 
and variance 

^ (3) 

where the overbar denotes an average. In this scenario 
the steps are balanced in the logarithmic space (the ex¬ 
pectation value of log{N) is kept fixed) but the arith¬ 
metic mean of N is growing in time, since exp{^) > 1. 
The reason for this is the asymmetry between growth 
and decline: since the per capita growth/decay rate is 
kept fixed for some time, the overall demographic benefit 
for the population includes not only the birth originated 
from the individuals that were present at s, but also from 
the individuals that were born between s and s -I- 1, and 
the opposite is true in the case of decline. The response of 
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the population to varying environmental conditions has 
an “inertia” that increases the overall demographic ben¬ 
efit during good times relative to the loss suffered during 
bad times, hence producing a net bias towards growth. 

However, in many scenarios the per-capita growth rate 
decreases when the population increases. For example, 
in a fixed size community with a zero sum game like the 
one considered in [26] , the fitness of a species determines 
its chance to replace an individual of another species by 
its own offspring; the more abundant a species is, the 
greater the chance of intraspecific competition so the in¬ 
ertia is weaker. In its extreme limit one can model this 
kind of behavior by taking Ng+i = "fsNs-, where js is 
again a balanced noise. Now it is the expectation value 
of N (rather than log{N)) which is kept fixed. We define 
these two types of environmental noise as logarithmically 
balanced (case A, where the ’’opposite” of iV —>■ N/2 is 
N —>■ 2N) and arithmetically balanced (case B, the ’’op¬ 
posite” of A —>■ N/2 is A —>■ 3A/2). Realistic systems 
with environmental stochasticity are, most likely, some¬ 
where between these two extremes, but the solutions we 
present below for these two scenarios provide the basic 
insights needed for consideration of the generic case, as 
explained in the discussion section. 


AGE-ABUNDANCE RELATIONSHIPS AND 
LIFETIME STATISTICS 


Within the above theoretical framework, we have de¬ 
veloped and solved a model of population dynamics un¬ 
der both environmental and demographic noise, making 
a distinction between logarithmically balanced (case A) 
and arithmetically balanced (case B) stochasticity. The 
results have been derived for a model with a geometric 
distribution of offspring, but, as explained in Appendix 
C, they are valid for any realization of the demographic 
noise provided that one is interested in the long-term be¬ 
havior of the system. 


Survival probability 

First we consider the survival probability, assuming a 
single individual (point speciation) at t = 0. It turns out 
(see Appendices B and C) that the survival probability 
at long times is governed by the equations: 

+ DS{t) - case A 

case B (4) 

In the limit 'y = D = 0 both equations reduce to Q, 
the equation obtained previously for pure demographic 
noise, as expected. Moreover, substituting 6 = 2DS/a^, 


changing the time units with t = t/D and rescaling the 
noise term to unit white noise, one finds a D indepen¬ 
dent equation. This implies that a plot of St ~ St vs. 
Dt = T should be independent of the strength of the 
environmental noise. 


From Eq. (161 one may extract the long-time asymp¬ 
totics of the survival probability. When the environmen¬ 
tal stochasticity is relatively weak the short time chance 
of extinction is controlled by the demographic noise. This 
leads us to propose that, for both the logarithmic and the 
arithmetic case, 


m 


T{tD) 

t 


(5) 


where the form of F{tD) is obtained analytically in Ap¬ 
pendix B II. 

In case A (log balanced) F[x) approaches yfxjr when 
X ^ oo and unity when a: —>■ 0. Accordingly, the chance 
of extinction at the t-th generation crosses over from the 
Galton-Watson universal limit (characterizing a process 
with pure demographic noise) t~^ to the first passage 
time asymptotics at t* ~ I/A. 

The intuition behind this result is clear. In the absence 
of demographic noise, the population preforms an unbi¬ 
ased random walk in the logarithmic space, hence it will 
survive until it reaches a threshold at, say, log{N) = 0 (a 
single individual). The theory of first-passage time for a 
Id random walker tells us that at long time the chance 
to survive decays like t~^/‘^. Since populations that stay 
alive for a long time typically reach high abundance, the 
long time behavior is controlled by this term while 
at shorter timescales environmental stochasticity is too 
weak and the behavior is controlled by the t~^ term of 
the purely demographic process. 

Figure]^ shows simulation results for the survival prob¬ 
ability in case A, i.e., the chance that the system did 
not go extinct until t, S(t). Indeed, this quantity decays 
like when the process is purely demographic, and it 
switches to t~^/'^ behavior at long times when the sys¬ 
tem is subject to environmental noise. Moreover, when 
t5{t) is plotted against Dt the data collapses as predicted 
above. 

In case B (arithmetic balance) F{x) approaches 
exp{—x/4:) when a; —>■ oo and unity when a; —)■ 0, as 
demonstrated in Figure]^ Again the short-term behav¬ 
ior is controlled by demographic noise, but the long-term 
survival probability is no longer a power law. The reason 
is that, unlike the log-balanced noise, in case B most of 
the species are shrinking in time (as can be seen easily 
by tracing the I/2-3/2 and 1/2-2 processes to the next 
generations). 

Returning to the empirical species lifetime problem, for 
a neutral theory with pure demographic noise the number 
of generations needed for a population to shrink to zero 
is typically its abundance, yielding unrealistic lifetimes 
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FIG. 2: A) A plot of St, the chance to survive at time t for a certain levels of overall diffusion coefficient D (log-log scale). The 
simulated process is the standard Wright - Fisher process with non - overlapping generations. B) Plotting St ■ t versus time 
in linear-log scale versus time; C) A plot of St ■ t against t ■ D. In this case the data collapses as predicted by Eq. 
collapse suggests that the proper scaling of the cumulative survival probability is St = with J-{x) ~ for 

J-{x) ~ const, when a; —>■ 0 


This 
a: —>■ 00 and 



FIG. 3: The survival probability of a species, 5, vs. time in 
case B (arithmetically balanced noise). The logarithm of tS{t) 
is plotted against t (inset), showing an exponential decay with 
D dependent slope. When plotted against Dt (main panel) 
the data collapse, indicating the dependence of J- on Dt. 


for common species. When environmental noise is intro¬ 
duced, the long-term dynamics (where long is relative to 
the strength of the environmental noise, i.e., t ^ ^/D) 
is dominated by cc = Dt, so a decent amount of envi¬ 
ronmental stochasticity will shorten significantly the life¬ 
time and will solve the problem with 0{N) generations 
scaling. Regarding the statistics of the species lifetime 
distribution, all the possibilities considered here - scal¬ 
ing with in case A, for pure demographic noise 
and exponential decay - have been suggested in the lit¬ 
erature based on the analysis of fossil data, see [22] and 
references therein. 


Age-abundance relationships 

Now let us turn to the dependence of abundance on the 
species’ age. In Appendix B we define and calculate Q{t), 
the mean abundance at time t of the ancestry originated 
from a a single individual at t = 0. However, when ex¬ 
amining an empirical community one considers only the 
species that have not yet gone extinct, so the relevant 
quantity for the age-abundance relationships is the aver¬ 
age size of a species at t conditioned on non-extinction, 
Nsurv(,t) = C{t)/S{t). As shown in Appendix B-III, there 
are strong differences between logarithmic and arithmetic 
noise, although the bottom line is similar. 

In case A (log balanced) C(0 ~ — l). Since 

iV.„™(t) ^ (6) 

Thus the abundance of the surviving species, instead of 
growing linearly, is growing exponentially in time. In Ap¬ 
pendix B-IV we explain that the average result presented 
here and the typical result differ from each other due to 
the skewness of the distribution, so in the typieal case 
the growth of Ngurv is subexponential 

NllUt) - (7) 

For case B, on the other hand, C = 1) since the noise 
in balanced in the real space. On the other hand 5 is 
decaying exponentially, so the net result is again, 

Nsurvit) ~ e^^/\ ( 8 ) 

DISCUSSION AND CONCLUSIONS 

The failure of the UNTB to account for dynamic pat¬ 
terns of populations and communities is known for a long 
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time m, and was stressed recently by many authors. Ba¬ 
sically, the OVn scaling of the demographic noise makes 
it inadequate to account for the observed fluctuations 
on all timescales. The tempo of the dynamics may be 
accelerated if one assumes a very large value of tr^ (as 
suggested, essentially, in m , see m) or by keeping the 
generation time as a free parameter (see, e.g., [21]), but 
any of these approaches carries its own difficulties. 

On the other hand, environmental stochasticity is 
known to be ubiquitous in living systems, affecting com¬ 
munities even under the most stable conditions (see, e.g., 
[29|b Incorporating this mechanism into the neutral the¬ 
ory is a required step in any case m- We have already 
showed that environmental noise increase substantially 
the heterogeneity of the species abundance distribution 
(see [55], Eq. 3), a feature that may account for the em¬ 
pirical results analyzed in [7] . The fact that this noise in¬ 
creases temporal fluctuations and decreases the timescale 
makes this project even more attractive. 

Under environmental stochasticity, the demographic 
rates of all individuals belonging to a species (roughly 
speaking, their fitness) are fluctuating coherently in time, 
and the species abundance varies accordingly. As ex¬ 
plained, the environmental noise may be ’’neutral” in two 
different senses. One scenario is when the relative fitness, 
when averaged over time, will be zero, this corresponds 
to logarithmically balanced noise or case A considered 
above. The other scenario, case B, occurs when the time 
average of the demographic gain is zero. 

As shown above, these two cases correspond to two dif¬ 
ferent species lifetime statistics. The chance of a species 
to survive decays like in case A and exponen¬ 

tially in case B. All these possibilities, —3/2 law, —2 law 
and exponential, were suggested for the tails of species 
lifetime distributions as extracted from fossil data [55] . 
More important is the transition of the general scaling 
from t (measured in units of generations) to Dt, allowing 
the environmental noise to control the extinction times. 

For the average abundance of a species at time t after 
point speciation conditioned on non-extinction, N^urv (t ), 
we obtained exponential growth (in case B) and a typi¬ 
cal stretched exponential growth (in case A) . This ap¬ 
pears to solve the “age-size” problem while the time 
to the most recent common ancestor scales with abun¬ 
dance in a purely demographic neutral model, it scales 
with the logarithm (or logarithm squared) of N in the 
presence of demographic stochasticity. Eqs. 0i suggest 
that the time from speciation to abundance N scales like 
log{N)/D generations (case B) or with log^ N/{2D) (case 
A). If, for example, one considers a set of 10® conspecihc 
trees for a frequent species in the Amazon basin (this 
is close to the contemporary figure, see recent survey in 
[31jb with about a 50y generation time, the neutral the¬ 
ory suggests 50 • 10® generations, more than the age of 
the universe. On the other hand, the left panel of Fig. 
[^suggests that,measured in units of a single generation 



FIG. 4: The mean of the logarithmic (red) and arithmetic 
(black) abundance fluctuations as a function of the abun¬ 
dance, extracted from simulation of a neutral community 
model with both demographic and environmental noise. A 
Moran-type continuous time process with discrete individuals 
was simulated. At each small time step {dt = 0.1) we gen¬ 
erate the number of births and deaths for each species. The 
number of births is a Poisson distributed number with mean 
ai{t)dt, where Qi(t) is the instantaneous growth rate (fitness) 
of species i. The number of deaths in the species is binomi- 
ally distributed, with the probability Ntot/K per individual, 
where Atot is the instantaneous size of the community and 
K — 10®. In addition, a Poisson number of new immigrants is 
drawn, with mean pdf, where /i — 2. Each immigrant founds 
a new species. The birth rates ai are given by ai = exp{'yi) 
where the 7 i are generated by an Orenstein-Uhlenbeck pro¬ 
cess, 7 i = — 7 i/T + 7 , so that they are Gaussian distributed 
with mean 0 and variance 0.001 and correlation time 2. This 
way the births rates are guaranteed positive. 

(50y) 13 ~ 0.1, so the neutral theory with environmental 
noise modifies the speciation time to about 50,000 years 
before present (case B) or about 150,000 ybp (case A). 
This, of course, is an extrapolation since the environ¬ 
mental noise may be either much smaller (if some of the 
short-term fluctuations are averaged out due to balanc¬ 
ing effects like the increase of species specific parasites) 
or much larger (if extreme events does not appear in the 
half-generation window considered), but the scales are 
clearly small enough to solve the age-abundance problem 
even if the estimation for D is modified by an order of 
magnitude. 

Finally, we would like to comment about the generic 
case. When a community is subject to a stabilizing mech¬ 
anism that keep its size fixed or almost fixed, the dynam¬ 
ics of a single species must be somewhere between the two 
extremes considered along this paper. Pure case A dy¬ 
namics cannot hold as there average population increases 
exponentially, while pure B dynamics is rejected since all 
species go extinct with probability one so the size of the 
community must shrink eventually to zero. In Figure [^ 
we plotted the mean of {Nt+At —Nt)/Nt (positive in pure 
case A, vanishes in case B) and the mean ln{Nt+At/Nt) 
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(vanishes in pure case A, negative in case B) vs. N as 
extracted from simulation of a community model with 
environmental noise (species specific randomly varying 
fitness). The situation is clearly between case A and case 
B - the arithmetic mean is slightly positive, the log mean 
is slightly negative - and shows some nontrivial N de¬ 
pendence. Our approach show how a rich behavior may 
arise in a simple model where we do not model interac¬ 
tion among different species. Further works will investi¬ 
gate the effect of environmental noise when considering a 
whole ecological community where species explicitly in¬ 
teract among each other. 
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APPENDIX A: SURVIVAL PROBABILITY IN A MODEL WITH PURE DEMOGRAPHIC NOISE 


We start by considering an event of point speciation, i.e., a new taxon that appears as one individual, the founder. 
Let us denote the chance that a single individual has n descendants in the sth generation hy Pn ■ In a constant 
population, perfectly neutral scenario with pure demographic noise, the number of offspring, n, that every given 
individual produces is independent of time and identically distributed, given by . The average number of offspring 
is unity, i.e., '^^nPn'’ = 1 and the strength of the demographic noise is characterized by ~ \)Pn'’ ■ 

This scenario was first analyzed by Galton and Watson (GW), who showed that such a population goes extinct with 
probability one [5S]. 

To generalize the GW result, one defines the generating function 


and its generalization to s generations. For example, P^'^ is the chance to have n grandchildren, while is the 

corresponding generating function. The successive generating functions are given by 

G(")(a:) = (9) 

By definition, = Pq'’ is the probability that the lineage originated from the founder went extinct by the sth 

generation, and G*^®“''^^(0) — G^^^(O) « dG^^\0)/ds determines the chance of extinction at the s generation. After 
many generations the chance of extinction is almost one, so G^®^(0) = 1 — ^s, ^ 1- Eq. § implies that 

1 - ^.+1 = ^ (1 - SsT « E (1 - , 


which leads to the recursive equation 


(5,+i =5s- (j'^51/2. 

Taking the continuum limit we then obtain that the survival probability is determined by the differential equation 

5{t) - 

This result implies that the long-time decay of the survival probability goes like S{t) « 2f{a‘^t), so the statistics of 
species lifetime admits a 1/t^ tail |55j. 


APPENDIX B: SURVIVAL PROBABILITY AND AGE-ABUNDANCE RELATIONSHIPS IN A MODEL 
WITH BOTH DEMOGRAPHIC AND ENVIRONMENTAL NOISE 

I. Neutral dynamics with environmental noise: The geometric neutral process 

In this section we derive our main results, implementing a model with a geometric distribution of offspring. As 
explained in the main text, two cases are considered here: a logarithmically balanced noise and an arithmetically 
balanced noise. 

In a neutral model with a geometric distribution but without any stochasticity, the chance of an individual to 
produce n offspring is P„ = 1/2"+^ and the corresponding generating function for a single generation is 

GW(x)=^a.-P„ = ^. 

n 

What makes this model easy to handle is the convenient structure of the generating function. The generating function 
for the population after s generations is obtained by iterating G^^\ and in this specific case the answer is immediate 

m. 


G^^\x) ^ = G(G(G...s times (G(a:)))) = " ^ 

n—0 ^ ' 


( 10 ) 
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Thus the distribution, apart from Pq^ remains geometric. From this one can easily derive the results given above for 
a purely demographic model for this special case of geometrically distributed births. 

To construct a model combining both demographic and environmental noise, we consider a discrete time dynamics, 
where for convenience we choose the time step to equal the generation time, so that all individuals of the sth generation 
reproduce and then simultaneously pass from the scene. 

Now let us consider the two random processes defined in the main text. The probability for each individual at the 
s-th generation to produce n offspring conditioned on an environmental noise determined 7 ^ is 

e'Ys'n 

= (e7, + l)r^+l ^ 


where case A is the logarithmically balanced noise and case B corresponds to arithmetically balanced noise. 

When 7 ^ = 0 for any s one obtains, of course, a purely demographic process. For nonzero 7 s the fitness (or the 
deterministic growth rate) of the population fluctuates, and D = 7^/2 characterizes the strength of the environmental 
noise. The model is neutral in the sense that 7 s is distributed identically for all species and so considered over long 
time scales, all species are demographically equivalent. 

The generating functions in the two cases are then. 


G(i)(x|7s) 
G(i)(a;|7s) = 


1 

" 1 + eT'= - xe-y- ’ 

1 

2 + 7 — a;(l + 7 ) ’ 


case A; 
case B. 


( 12 ) 


Using the recurrence relation Eq. Q and the offspring generating function, Eq. (12), one can obtain [32] the 
general form of G(®^(a;| 7 s) as 


gW(®|7s) = 


dg “)“ bg • X 


Cs +da ■ X’ 

where the values of the constants aa,bg,Ca and dg satisfy [for a log-balanced noise (case A)] the recurrence relation 


(13) 


ffls + l 


0 

1 

0 

0 


dg 

Cs + 1 


_g7a + l 

g7s+l _|_ 1 

0 

0 


Cg 

&S + 1 


0 

0 

0 

1 


bg 

dg+1 


0 

0 

_g73+l 

g7s + l 1 


ds 


(14) 


with the initial conditions ai = 1, bi — 0, ci = -I-1 and di = —The corresponding equations in case B are 

obtained by replacing e'’' by 1 -I- 7 in all the above expressions. 

The chance of survival until the sth generation ds = 1 — G'’(0) is simply 


Sg — 1 dg^Cg. 


Accordingly, the survival probability satisfies the stochastic recursion relations: 


e-^» -f Ss ’ 

^ (1 + 7)<^^ 

1 -I- (1 -I-7)(5s’ 


case A; 
case B. 


(15) 


When 5s ^ 1, it satisfies 5s+i = in case A and 5s+i = (1 + 7 )'^^ ~ <^5 in case B. One has to be careful in 

translating the equation for case A to a stochastic differential equation, as exp{'js) has a non-zero expectation value, 
namely cosh 'J2D. Thus, in terms of 5 there is an extra term driving 5 to larger values. Taking this into account, 
in the weak environmental noise limit, 7 s 1 , one gets the stochastic differential equations that were quoted in the 
main text: 


dS{t) 

dt 




case A; 


dt 


case B. 


(16) 
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FIG. 5: The Morse potential V{y) = exp(2y)/4D — exp(y)/2 for D = 0.1. 


II. Long-time asymptotics 


The Langevin equations (16 1 should be interpreted in the Ito sense, since they were derived as the white noise, 
continuous time limit of a non-overlapping generation model that assumes zero relaxation time, i.e., that the population 
follows the instantaneous growth rate determined by 7 [331 HI] ■ When the relaxation time of the population or the 
community is finite (i.e., when the demographic rates respond slowly to the changing environment, with respect to the 
noise correlation time) the effective strength of the environmental noise amplitude decreases, but the interpretation 
is still Ito. Accordingly, in the logarithmic space y = ln{5), Eqs. (16) takes the form: 


ds 


dy , , , , 

— = 7 (s) - exp(y), 

= 7(s) - D- exp(?/). 


case A; 
case B. 


(17) 


It is clear, now, that as i5 —>■ 0, i.e., y —> — 00 , in case A the system performs an unbiased random walk in the log 
space, hence in the long run one expects that the probability of a taxon to have a lifetime t will behave like i.e., 

that 5{t) ^ In case B, on the other hand, the random walk is biased to the left, and the survival probability 

should decrease exponentially in time. 

To demonstrate that, we notice that the corresponding Fokker-Planck Equations are 


dsP{y, s) = DdyP{y, s) + dy (exp(y)P(?/, s)) , case A; 

dsP{y, s) = DdyP{y, s) + dy {[D + exp{y)]P{y, s)) , case B. 


(18) 


With the substitutions P{y,s) = exp{—e^/2D)'tjj{y, s) (A) and P{y,s) = exp(—j//2 — e'^/2D)ip{y, s) (B) one gets 


tp{y, s) = Dtp'' + 


/exp(y) exp( 2 y)\ 




4D 


- 


V', 


tp{y, s) = Dtp" - 


\ AD ^ 4 ' 


case A; 
case B. 


(19) 


These are Schrodinger equations in imaginary time with an exponentially decaying potential (case B) and a Morse 
potential (case A, see Fig. [^. 

A wavepacket which is initially localized at small negative values of y (corresponding to 5 ~ 1) will eventually reach 
the region of large negative j/ (<5 —>■ 0 ), where the potential is negligible and the motion is almost purely diffusive. 
This is clearer in case B, where the long-term behavior of tp (neglecting the exponential term) is. 




e iDt e 


t 


( 20 ) 


so the wavefunction diffuses in the log space, but this diffusion is superimposed of an exponential decay exp(—i7t/4) 
(we have replaced the generation parameter s by time t, as in the long-term the changes of 5 over a single generation 
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FIG. 6: Numerical simulation of the Fokker-Planck Eq. 
Dt = 10 and Dt = 25. 


(161, compared to the analytic approximation, Eq. (211, for D = 0.1, 


are small). For case A an asymptotic matching approach leads to the uniform approximation 


P{S,t) 


^ ^-(ln(SID)+'yE)ViDt^-S/D 

S^/ irDt 


( 21 ) 


Where je is the Euler constant. This approximate solution is shown in Fig. together with a direct numerical 
solution of the Fokker-Planck Eq. We see that the agreement is very good, and improves with time. At long 

times, the expectation value of i5, 5 « y/Djn, so that the extinction rate has a tail. For case B we have. 


P{6,t) 


1 (in6 + D) 

_g 4Dt 

5y/'KDt 


( 22 ) 


so that the log-normal distribution moves toward negative values of ln[5) at a constant rate, leading to a typical <5 
which decays exponentially in time. 


Ill: Abundance-age relationship 


We now turn to analyze the effect of environmental noise on the abundance-age relationship in neutral theory. As 
explained in the main text, there are two ingredients that determine the age-abundance relationships: the chance of 
a species to reach abundance N after s generations, and its chance of survival. The ratio between these factors gives 
Ngurv, the average size of a surviving species. The chance of survival was calculated above, and now we will calculate 
the average abundance. 

In case B the result is quite trivial: since the demographic and the environmental stochasticity are both balanced 
in the abundance space, the average abundance of a species is fixed. Accordingly, we perform the analysis here for 
case A, and will e xplain at the end how the results relate to our prediction for case B. 

From Eq. (13), we can see that the width of the geometric distribution of the abundance at time s is given by 


^(s) = —1/ln(—ds/cs). Note that ds and Cs satisfy the same recurrence relations, and the only difference between them 
is the initial condi tion s. This allows us to consider only dg and then to carry over the results to Cs and hence calculate 
their ratio. From (141 it is clear that ds satisfies 


The solution to this is easily verified to be 


Given di = — "w® can write this as 


e^°(ds - ds-i). 

(23) 

s k 



(24) 

fc=l j=2 


s k 


EH'” 

fe=ii=i 

(25) 


As explained above, the equations for c is the same, except that ci = 1 -I- e'^'b so that 


s k 

k=0j=l 


( 26 ) 
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FIG. 7: A gr^ 
function Eq. (141 


of the ensemble-averaged width ^{s) calculated by averaging as computed 
over 10® realizations of binary noise 7 = ± 70 , with 70 =0.1. Also show in the 


from the analytic gen erat ing 
analytic formula, Eq. (27l. 


Thus, ds = —(cg — 1) and so 1/C = — ln(—dg/cg) = — ln(l — l/cg). It should first be noted that for s 1, U <C 1, Cg 
is large. For the case of binary noise, for example, 7 = ± 70 , the minimum possible value of Cg is for the case where 
7 g is always negative, in which case Cg « l/ 7 o- Thus, C ~ Cs, so that 


k—Oj=l 


E' 

k=0 


,Dk/2 _ 


D 


(^^Ds/2 _ 


(27) 


So we see that Cg is exponentially large for large s. 

For case B, as explained, one should replace e'^ by 1 -I- 7 , and the trivial result is a linear growth, C ^ 


t. 


The result (27) has two limits: when s 1/71, the environmental noise is negligible and the average abundance 
of a species (conditioned on its presence in the community) grows linearly in time. However, at long times the 
environmental noise controls the system and the typical size of surviving species grows exponentially with s. In this 
regime the typical timescale (in generations) needed for a species to reach abundance N will scale like log{N)/D, as 
opposed to the linear N scaling for pure demographic noise. 


IV: Average vs. typical age-abundance relationships in case A 


Our result for case A, Eq. (27), might appear at first glance to be in contradiction with the results proved by Geiger 
and collaborators 


]. These authors show that, in case A 




(28) 


and decided that this log-balanced scenario is “critical”, as opposed to case B (or, in general, any noise which gives 
a negative bias in the log-space) which is subcritical. 

To explain this apparent contradiction, let us stick to the method used in [35], averaging first over demographic 
noise. In the log-abundance space, what one has after taking this average is a random walker that starts at zero 
(one individual) and moves randomly with no bias until s. Clearly, the population at s, conditioned on survival 
and averaged over demographic noise, is exp{y), where y is the location of a random walker that starts at s = 0 
and subsequently never crosses zero, since otherwise it goes extinct. Such a constrained random walk is known as a 
meander. To calculate the probability that the meander is at y at time s one solves for a RW that start at Xo using 
images, taking xq to zero at the end. The result is 


y^-y^/ADs 

5372 


P(j/,s) = A 


(29) 
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Accordingly, the maximum likelihood for the position of the meander at s is (for large s), 

Yml = V2Ds. 


(30) 


This implies that the typical size of a population at s, conditioned on non-extinction is, 

NlyP, « e*' = exp(V^). 

On the other hand, the average value of will be, 

= -y^/4Ds 




dy 


g3/2 


(31) 


(32) 


Note that the result for the average comes from the peak of the integrand at y* = 2Ds (this is a Laplace integral) and 
P{y*) is exponentially small as s —>■ oo, so this result is consistent with [35j since the average comes from exponentially 
rare events (with probability one this will not be the case for any specific history). On the other hand, P{yML) is 
©(l/t), which is not negligible. However, the typical growth is subexponential, and clearly agrees with (p^. 


APPENDIX C: UNIVERSALITY 


Throughout section B we have studied the geometric neutral process because of its convenient properties. Of course, 
any result that depends on the specific properties of a certain distribution cannot be relevant to the generic case, 
where the distribution of number of descendents per individual is, in most cases, unknown and there is no reason to 
believe that it belongs to any particular simple distribution. 

To illustrate the generality of our results, let us consider the survival probability dg. In the generic case one may 
define a probability distribution function P„(s), the chance of an individual to produce n offspring during the s 
generation, and of course the two important summary statistics that characterize this distribution are its mean (in 
case A, h = expijs), in case B, h = 7 s and its variance. Here we assume that the environmental noise is weak so 
one can neglect the variance fluctuations and varin) = cr^, where tr^ is the variance of the purely demographic model 
(cr^ = 1 for Poisson distribution with average one, = 2 for the geometric distribution and so on). 

The generating function recursion relation implies that 




In the long time limit the chance to survive is small, so G^^\x) = 1 — ^s- Plugging this into (33) one gets: 


1 - (5s 


“ <5(5-i)” « + 2 ■ 

n T) ' 


(33) 


(34) 


Accordingly: 

2 

(5 « ^s — ^s-i = (e^” — 1)(5 — case A 

2 

6 = jsd — ^<5^ case B. (35) 

The second equation is what appears in Eq. 4 of the main text. In case A one needs to expand the exponent in 7 
and to replace 7^/2 by D to obtain the equation presented in the main text. 







